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.^ I The interplay rich between algebraic geometry and string and gauge the- 

ories has recently been immensely aided by advances in computational al- 

(3 ■ gebra. However, these symbolic (Grobner) methods are severely limited by 

Cn I algorithmic issues such as exponential space complexity and being highly 

sequential. In this paper, we introduce a novel paradigm of numerical 
algebraic geometry which in a plethora of situations overcomes these short- 

^ I comings. Its so-called 'embarrassing parallelizability' allows us to solve 

many problems and extract physical information which elude the symbolic 
methods. We describe the method and then use it to solve various problems 
arising from physics which could not be otherwise solved. 
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1 Introduction 

Increasingly is the field of computational algebraic geometry becoming a signif- 
icant tool to contemporary theoretical physics [1]. This is becoming ever more 
transparent especially in string theory and in gauge theory; indeed, the plethora 
of vacuum configurations in string and M-theory is clearly best suited for large- 
scale intensive computerization 0, so too can new hidden structures in gauge 
theory be uncovered by computing new geometric invariants (cf. e.g. [3H5]). 

The rapid progress in algorithmic geometry as well as computer power has 
been a crucial ingredient to such advances. Key to much of the algorithms an- 
alyzing simultaneous multi-variate polynomial systems of high degree, to which 
the aforementioned physical situations can ultimately be reduced, is the Grobner 
basis technique (for an excellent review in this context, cf. |6]). Roughly speaking, 
given the ideal generated by the vanishing set of the polynomials, the so-called 
Buchberger Algorithm (BA) or its refined variants can compute a new system of 
equations, called a Grobner basis [7]. 

The BA reduces to Gaussian elimination in the case of linear equations. Sim- 
ilarly it is also a generalization of the Euclidean algorithm for the computation 
of the Greatest Common Divisors of a univariate polynomial. Recently, more 
efficient variants of the BA have been developed to obtain a Grobner basis, e.g., 
F4 [8], F5 [^ and Involution Algorithms jlO]. Symbolic computation packages 
such as Mathematica, Maple, Reduce, etc., have built-in commands to calculate 
a Grobner basis. Moreover, Singular |TT], COCOA jT^ and Macaulay2 [13] are 
specialized packages for Grobner basis and Computational Algebraic Geometry, 
available as freeware, and MAGMA |14] is also such a specialized package available 
commercially. 

In [S1[I51[I5], the power of Grobner basis methods was harnessed for the sake 
of answering questions pertinent to string and particle phenomenology. A pub- 
licly available computational package, called Stringvacua, is designed to interface 
Mathematica with Singular and has built-in utilities for string phenomenology. By 
extracting important information such as the dimension of the ideal, the number 
of real roots in the system, etc., many issues such as stability or supersymmetry 
of the potential or the branches of moduli space of vacua, can be settled using 
only a regular desktop machine. 

However, even with such reassuring techniques, there are a few problems with 
symbolic methods: the BA is known to suffer from exponential complexity, both 
the RAM required by the machine and the computation time increases exponen- 
tially with the number of variables, equations, degree and terms in each polyno- 
mial. Thus, often for even seemingly small sized systems, one may not be able 
to compute a Grobner basis. Furthermore, BA is also usually less efficient for 
systems with irrational coefficients and one habitually has to resort to random- 
izing over the space of integer (or prime) coefficients and work over finite fields. 
Another drawback is the highly sequential nature of the BA and thus efficient 



parallelizations currently do not exist. 

When the information we are after is not sophisticated quantities such as coho- 
mology of vector bundles, which, for example is central to understanding particle 
spectra and Yukawa couplings in heterotic compactifications, it is natural to won- 
der whether one could appeal to numerics. Indeed, in a typical situation when 
one is confronted with the minimization of a potential to see how many isolated 
solutions there are, what are their approximate values and whether supersym- 
metry is preserved, numerical recipes specifically tailored to polynomial systems 
would seem ideal. 

The purpose of this paper is to explain precisely such a method, called the 
Numerical Algebraic Geometry (NAG) and to introduce it in our context. 
The method was first used in lattice field theory [TTj where all the stationary 
points of a multivariate function called the lattice Landau gauge fixing functional 
|18f[2T] needed to be found. We will see how the shortcomings of the Grobner 
basis methods may be overcome with this numerical method. 

In general, the problems we encounter fall into two categories: (1) there is a 
finite number of isolated extrema and (2) there is a continuum of solutions. The 
first constitutes zero- dimensional ideals and arises, for example, when finding the 
number and numerical positions of minima of a given potential, be it from a het- 
erotic string model or from a cosmological background. The second gives us an 
afiine algebraic variety and involves queries from an algebraic-geometric perspec- 
tive. These could come from analysis of the space of vacua of scalar fields in a su- 
persymmetric gauge theory, or from complex geometries arising from AdS/CFT. 
The central method in NAG is the numerical polynomial homotopy continuation 
(NPHC): it can be shown that for a given system of polynomial equations to be 
solved, a homotopy between the system and a new system (which is easier to solve 
and share many features with the former system) can be constructed (see [2^123] 
for a detailed description of this method in this context). Then, one tracks paths 
starting from each solution of the new system as one moves towards the original 
system along the homotopy, to finally obtain all the solutions of the original sys- 
tem. The NPHC method is used to find all the complex (which obviously include 
real) solutions of a system which is known to have only isolated solutions, i.e., the 
first category of the systems. The NPHC method has been immensely useful in 
various areas of theoretical physics including statistical mechanics, particle phe- 
nomenology, string phenomenology, lattice field theories, etc. [T7 t fT8 |[22] - [28] . For 
the second category of systems, the NPHC method can be extended to describe 
the solution space completely (i.e., dimension, degree, and number of irreducible 
components). 

Indeed, one of the greatest computational hurdles to Grobner techniques is 
the primary and irreducible decomposition of ideals which gives crucial data of 
the branches of the moduli space. We will address ample examples from both 
categories, demonstrating how many questions can be settled by NAG. 

The organization of the paper is as follows. We begin in § |2] by describing 



the NPHC method, first pedagogically for a univariate example before gener- 
alizing it to the multivariate polynomial ideal. We then apply the method to 
the two categories of circumstances. In § |3l we treat several systems from het- 
erotic and M-theory models which the symbolic methods are already known to 
be prohibitively difficult. Next, in § |U we discuss how to perform numerical 
irreducible decomposition to systems which are positive dimensional, and apply 
to four different situations arising from string and gauge theories. We conclude 
with prospects in § O 

2 The Numerical Polynomial Homotopy Contin- 
uation Method 

In this section, we will outline one of the most powerful tools in numerical 
algebraic geometry, the so-called numerical polynomial homotopy continuation 
method, or NPHC. Since this method is already introduced in particle physics 
and statistical mechanics areas in |17 t [TS | 1221 - 125) . we will only give a brief ex- 
planation as to the principles and algorithms involved as well as advertise some 
available softwares. We note that in these references only one kind of homotopy, 
called the Total Degree Homotopy, was used whereas in the present paper we 
introduce two new kinds of homotopies, the 2-Homogeneous Homotopy and the 
Polyhedral Homotopy, which in many ways are more powerful than the former. 

We know from the Fundamental Theorem of Algebra that for a single variable 
polynomial equation of degree k G Z>o, say f{x) = Yli=o^i^\ with coefficients 
Oj and the variable x both defined over C, the number of solutions is exactly k if 
ttk 7^ 0, counting multiplicities. To solve a univariate equation, many numerical 
methods such as the companion matrix trick for low degree polynomials, and 
divide-and-conquer techniques for high degree polynomials, are available. Here, 
we introduce the NPHC method, which can then be extended to the multivariate 
case in a straightforward manner [29|l30]. 

Briefly, in the NPHC method, flrst one writes down the equation or system 
of equations to be solved in a more general parametric form, then solves this 
system at a point in parameter space where its solutions can be easily found, 
before flnally tracking these solutions from this point in parameter space to the 
point in parameter space corresponding to the original system/problem. This 
approach can be applied even to non-algebraic equations which exhibit a contin- 
uous dependence of the solutions on the parameters. However, there exist many 
difficulties, which do not arise in the polynomial setting, in making this method 
a primary candidate method to solve a set of non-algebraic equations. 

Let us consider the univariate equation z"^ — 5 = 0, pretending that we do not 
know its solutions, i.e., z = ±-\/5- We thus begin by defining the more general 



parametric family 

H{z, t) = (1 - t){z^ - 5) + -it{z^ - 1) = 0, (2.1) 

where t G [0, 1] is a parameter and 7 is a complex number. For t = 1, we have 
z^ — 1 = and at t = we recover our original problem. The problem of getting all 
solutions of the original problem now reduces to tracking solutions of H[z, t) = 
from t = 1 where we know the solutions, i.e., z = ±1, to t = 0. The choice of 
z^ — 1 in Eq. (12.11) . called the start system, is made for the following reasons: this 
system has the same number of solutions as the original problem and is easy to 
solve. For multivariate systems, a clever choice of a start system is essential in 
reducing the computation as we shall see. 

We briefly mention the numerical methods used in path-tracking from t = 1 
to t = 0. One of the ways to track the paths is to solve a differential equation 
that is satisfled along all solution paths, say z*{t) for the i*'* solution path, 

dH{z*{t),t) _ dHiz*it),t)dz*{t) dH{zm,t) _ 

dt dz dt dt ■ ^ ' ' 

This equation is called the Davidenko differential equation. We can solve 
this initial value problem numerically - again, pretending that an exact solution 
is not known - with the initial conditions as 2*(0) = 1 and ^2(0) — ~^- (Since 
we also know that H{z*{t),t) = 0, predictor-corrector methods, such as Euler's 
predictor and Newton's corrector, are used in practice to solve this initial value 
problem.) We shall not discuss the actual path tracking algorithms further, but it 
is important to mention that in these algorithms are designed to handle almost all 
apparent difficulties such as tracking singular solutions, multiple roots, solutions 
at infinity, etc. 

We now return to the complex number 7, where we will consider 7 = e*^ with 
6 G [0, 27r) chosen generically. It is known that for all but finitely values of 6* G 
[0, 27r), the paths are well-behaved for t G (0, 1], i.e., for the whole path except the 
end-point [29|. This makes sure that there is no singularity or bifurcation along 
the paths. This is a remarkable technique, called the 7-trick in the literature, 
and constitutes the reason why NPHC is guaranteed to find all isolated solutions. 

There are several sophisticated numerical packages well-equipped with path 
trackers such as Bertini [31], PHCpack [32], PHoM [33], HOMPACK [31] and 
HOM4PS2 [35|l36] . They all are available as freewares from the respective research 
groups. In the above example, PHCpack with its default settings yields the solu- 
tions z = -2.23606797749979 + z 0.00000000000000 and z = 2.23606797749979 + 
z 0.00000000000000. Thus, it gives the expected two solutions of the system with 
a very high numerical accuracy. 

2.1 Multivariate Polynomial Homotopy Continuation 

We will now generalize the NPHC method explained above to the multivariate 

case, say P{x) = 0, where P{x) = {pi{x) , . . . , Pm{x)) and x = {xi, . . . ,Xm), 
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that is known to have isolated solutions (i.e., a O-diinensional ideal). Of course, 
in many cases where one does not know the dimensionahty of the solutions to 
start with, one has to check if the given system is an 0-dimensional first. In the 
following sections, where we explain the extension of the NPHC method to the 
positive dimensional systems, we will describe a concrete way of finding if the 
system possesses only isolated solutions. One quick and dirty way of seeing that 
a square system (i.e., number of equations and number of variables are the same) 
has only isolated nonsingular solutions is to add the determinant of the Jacobian 
matrix to the system and verifying that the combined system has no solutions. 
In the 0-dimensional case, the system has maximal rank at each solution. Having 
said this, in many physical applications, we know a priori that the given system 
has only isolated solutions. 

We can construct a homotopy, or a set of problems, similar to the aforemen- 
tioned one-dimensional case, as 

H{x,t) = {l-t)P{x)+-ftQ{x), (2.3) 

where Q{x) is a system of polynomial equations, Q{x) = {qi{x), . . . , qm{x)) with 
the following properties: 

1. The solutions of Q{x) = H{x, 1) = are known or can be easily obtained. 
Q{x) is called the start system and the solutions are called the start solu- 
tions. 

2. The number of solutions of Q{x) = H{x, 1) = is equal to an estimated 
number (or an upper bound) of the solutions for P{x) = 0. 

3. The solution set of H{x, t) = for < t < 1 consists of a finite number of 
smooth paths, called homotopy paths, each parameterized by t G (0, 1]. 

4. Every isolated solution of H{x, 0) = P{x) = can be reached by some path 
originating at a solution of H{x, 1) = Q{x) = 0. 

We can then track all the paths corresponding to each solution of Q{x) =0 
from t = 1 to t = and reach H{x, 0) = P{x) = 0. By implementing an 
efficient path tracker algorithm, we can get all the isolated solutions of a system 
of multivariate polynomials just as in the univariate case. There are several 
upper bounds on the number of solutions of the system P{x) = 0, which yield 
alternatives to constructing Q{x) and the homotopy H{x,t) for the multivariate 
case. The ones that we will consider are the (1) Total Degree Homotopy, (2) 
2-Homogeneous Homotopy, and (3) Polyhedral Homotopy. We introduce the 
Total Degree Homotopy below and discuss the 2-Homogeneous and Polyhedral 
Homotopies in Appendices. 



2.1.1 Total Degree Homotopy 

The Total Degree Homotopy arises from the Classical Bezout Theorem, which 
asserts that a system of m polynomial equations in m variables has at most 
Y\T=i ^i isolated solution in C", where di is the degree of the ith polynomial. This 
bound, called the Classical Bezout Bound (CBB), is sharp for generic values (i.e., 
roughly speaking, non-zero random values) of the coefficients. The genericity is 
well-defined and the interested reader is referred to |29| for details. The homotopy 
constructed using the CBB is called the Total Degree Homotopy. The start system 
Q{x) = can be taken for example as 

f xt-l\ 

= 0, (2.4) 



Q(x) 



xi^ -1 



\xt-l ) 



where di is the degree of the i*^ polynomial of the original system P{x) = 0. 
Eq. 02. 4p is easy to solve and guarantees that the total number of start solutions 
is YYiLi ^j) ^^^ of which are nonsingular. The Total Degree Homotopy is a very 
effective and popular homotopy which is used in actual path trackers. 

The advantages of the Total Degree Homotopy are (1) the CBB is easy to 
compute, and (2) the start system based on the CBB can be solved quickly. The 
drawback of it is that the CBB does not take the sparsity of the system into 
account: systems arising in practice have far fewer solutions than the CBB, so a 
large portion of the computational effort is wasted. 

The 2-Homogeneous Homotopy is constructed by first writing C"* = C'^ x 
(^k-m j?Qj, some < A; < ?Ti, which is accomplished by partitioning the original 
variables into two groups. This has the advantage of incorporating some of the 
structure of the given polynomial system P[x) into the start system Q{x). The 
corresponding bound, called the 2-Homogeneous Bezout Bound (2HomBB), is 
often tighter than the CBB when the polynomial system P{x) has a naturally 
arising partition of the variables, which occurs in the examples below. Given a 
partition, the 2HomBB is easy to compute and the start system can be solved 
quickly via linear algebra. 

The Polyhedral Homotopy uses the monomial structure of the given polyno- 
mial system P{x) based on the Bernstein-Khovanskii-Kushnirenko (BKK) The- 
orem [37H39] to yield the BKK bound. Essentially, this upper bound on the 
number of complex solutions is obtained by computing the mixed volume of the 
convex hull of the Newton polytope (which is based on the exponents of the 
monomials appearing) of each equation. Since we must introduce some jargon to 
fully describe the 2-Homogeneous and Polyhedral Homotopies, we do it in Ap- 
pendix |X] and [B], respectively. We note that, as with the CBB, the 2HomBB and 
BKK bound are also generically sharp with respect to the family of polynomial 
systems under consideration. 
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Since the BKK count utilizes sparsity of the system, the BKK bound is tighter 
than the CBB and 2HomBB resulting in fewer honiotopy paths to track. The 
apparent drawback is that to compute the BKK bound and solve the start system 
itself, the Polyhedral Homotopy requires additional computational effort. As we 
will see, the difference between the CBB and BKK bound suppress this drawback 
of the Polyhedral Homotopy by saving a large portion of the computational effort 
over the Total Degree Homotopy. The relationship between the 2HomBB and 
BKK count vary based on the system. 

2.1.2 Parallelizability of the NPHC Method 

While going through the algorithm for the NPHC method, one could have noticed 
that to track paths from t = 1 to t = 0, i.e., from the start system to the original 
system, each start solution does not need any knowledge about any other path. 
In other words, each solution path can be tracked completely independently of 
all others. This feature makes the method embarrassingly parallelizable. This 
is sharply different than the BA which is known to be a sequential algorithm, 
i.e., each step requires the knowledge of the previous step. Thus, the NPHC 
method is destined to be much more efficient to solve bigger systems which are 
way beyond the reach of the Grobner basis techniques. 

2.1.3 Numerical Solutions of Multivariate System 

While tracking all the paths using either Total Degree, 2-Homogeneous, or Poly- 
hedral Homotopies, only those paths leading to solutions of the original system 
will converge at t = 0, that is, some of the paths will diverge. The values of the 
variables at t = 0, i.e., x{t = 0), after the path-tracking are the solutions of the 
original system. Obviously, they will be in the numerical form. It is useful to 
define what we mean by a solution here. For the multivariate case, a solution is 
a set of numerical values of the variables which is within a given tolerance Ag^j 
(we will take ~ 10~^° in our ensuing calculations) of an actual solution. Since the 
variables are allowed to take complex values, all the solutions come with real and 
imaginary parts. A solution is a real solution if the imaginary part of each of the 
variables is less than or equal to a given tolerance, A]r (~ 10^^ is a robust toler- 
ance for the equations we will be dealing with in the next section, below which 
the number of real solutions does not change). All these solutions can be further 
refined to within an arbitrary precision (up to the memory and computational 
limits placed on by the machine). 

The obvious question at this stage would be if the number of real solutions 
depends on Ajg. To resolve this issue, we use a recently developed algorithm 
called "alphaCertified" which is based on the so-called Smale's a-theory. This 
algorithm certifies which nonsingular solutions of the given polynomial system 
are real using either exact rational arithmetic and arbitrary precision fioating 



point arithmetic |40| . This is a remarkable step, because using alphaCertified, we 
can prove that a solution classified as a real solution is actually a real solution 
independent of Ajr, and hence these numerical approximations are as good as the 
exact solutions. 

3 Isolated Solutions and Phenomenology 

We will now apply the principles introduced in the preceding section on NPHC 
to questions which arise in theoretical physics, especially in string phenomenol- 
ogy and supersymmetric gauge theory. Our discussions will roughly fall into two 
categories to which we alluded in the discussions above: (1) polynomials systems 
which have any isolated solutions; these are clearly of importance when dealing 
with issues such as finding and stabilizing vacua given some effective potential, 
and (2) systems which can have solution components; these are crucial to un- 
derstanding various geometrical properties of gauge theories such as the moduli 
space of vacua or Calabi-Yau spaces in the AdS/CFT context. We will first dis- 
cuss (1) and then turn to (2) in the next section, illustrating in each case with 
ample examples, particularly those which defy more standard techniques. 

Much of current research in string phenomenology is focused on developing 
methods to find and analyze vacua of four dimensional effective theories for su- 
pergravity descended from flux compactiflcations. Stated in explicit terms, one 
is interested in flnding all the vacua (usually, isolated stationary points) of the 
scalar potential V of such a theory. In particular, given a Kahler potential K, 
and a superpotential W, for uncharged moduli flelds, the scalar potential can be 
obtained as 

V = e^[K^^ DaW DbW-3\W\'^] , (3.5) 

where Da is the Kahler derivative (9^ + OaK and K^^ is the inverse of Kab = 
OaObK. Once the vacua are found, one can then classify them by either using 
the eigenvalues of the Hessian matrix of V or by introducing further constraints 
such as 1^ = 0. 

Finding all the stationary points of a given potential V amounts to solving 
the stationary equations, i.e. solving the system of equations consisting of the 
flrst derivatives of V with respect to all the flelds equated to zero. The stationary 
equations for V arising in the string phenomenological models are usually non- 
linear. In the perturbative limit, W usually has a polynomial form. This is an 
important observation since we can then use algebraic geometric concepts and 
methods to extract information about V. 

For the systems known to have only isolated solutions, we are in the situa- 
tion of 0-dimensional ideals, a Grobner basis using a lexicographic ordering of 
the monomials always has at least one univariate equation and the subsequent 
equations consist of an increasing number of variables, i.e., it is in a triangular 
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form. The solutions of a Grobner basis is always the same as the original sys- 
tem, but the former is easier to solve due to its triangular form as the univariate 
equation can be solved either analytically or numerically using straightforward 
methods. Then, by back-substituting the solutions in the subsequent equations 
and continually solving them, we can find all the solutions of the original system. 

3.1 Some Toy Examples 

One of the main features of the "Stringvacua" package [151 [IS] is to address 
problems with a 0-dimensional ideal using Grobner basis techniques. Let us 
first take some of the toy examples discussed therein and see how our numerical 
techniques can extract some of the requisite information without resorting to the 
often expensive Grobner algorithms. 

Sysl: A Single- Modulus Example: We begin with a single- modulus toy 
example. Let the Kahler potential K and superpotential W be given as 

K = -3 log(T + T) , W = a + bT^ . (3.6) 

Note that the field T comes along with its complex conjugate. So even though 
they can be treated as different variables by merely relabeling them, they are 
not actually independent variables. To avoid this problem, we can write them in 
terms of real and imaginary parts, i.e., T = t + ir where t and r are real. The 
potential is 

V = ^{4b{5b{t^ + T^y - ?,a{f - 21tV2 + 35tV^ - 7t^))) (3.7) 

Ob 

which has 2 variables. To find the stationary points of V , we need to compute 
the zero locus of the partial derivatives of V with respect to variables t and r: 

OT 7" -I 

^ = (46(56(13t2-r2)(t2 + 7-2)6 _3a(5t6-63tV2 + 35tV + 7r^)))=0, 

at 3t^ 

^ = -(566r(56(t2 + T^f + a(9t^ - 3Qt\^ + 9r^))) = . (3.8) 

UT OL 



For general values of parameters a and b it was already known j6l[15] that the 
system becomes prohibitively difficult to analyze using symbolic methods. For 
now, we take a = b = 1. 

We also note that the stationary equations in this example involve denomi- 
nators. Since we are not interested in the solutions for which the denominators 
are zero, we clear them out by multiplying them with the numerators appropri- 
ately. In these equations, all the denominators are multiples of t. The condition 
that none of the denominators is zero can be imposed algebraically by adding 
a constraint equation as 1 — 2;t = with z being an additional variable. Thus 
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there are now 3 equations in 3 variables. One could impose this by the method 
of saturation |15| . 

Stringvacua can easily solve this system: this system has 6 real solutions. Note 
that this computation uses the Grobner basis technique implemented in Singular. 

We now turn to solve the above system using Bertini and HOM4PS2. The 
CBB of this system is 364, while the 2HomBB, using {t, r} x {z}, and the BKK 
bound are both 182. In the end, there are 56 finite nonsingular complex solutions, 
of which six are real. Thus, we have produced the results known from symbolic 
methods by using the NPHC method. Note that the choice of a and b was crucial 
here. For a generic choice for a and b, the computation using symbolic methods 
via Stringvacua becomes extremely difficult already for this small system, whereas 
the NPHC method takes the same time to solve it as before. Since we have all 
the real solutions, we can then compute the Hessian of V at the real solutions 
and separate out the physically interesting vacua. 

Sys2: Two-Moduli Model: Consider the Kahler potential and superpoten- 
tial 

K = -3\og{T + f)-\og{S + S), 

W = aS + bST + cT\ (3.9) 

with two fields T = t + it and S = s + ia. Again, a, b and c are the parameters 
and for convenience we chose them to be 1, —1, 1, respectively. The potential is 

V = _L^(-5t4 + 3(3^ + ^2) -2tV + 3r^- (s2 + a2)(5t2-3r2) -6(5st2- 
t{s^ + CT^) + 2taT - st"^) + 2(13st^ + tVr + Ostr^ - 3crr^)) (3.10) 

and the stationary equations are 

= ^(-5t4-9(s2 + a2) + 2tV-9r^ + (s2 + a2)(5t2-9r2) 

+6(5st2 - 2t(s2 + a^) + Atar - ?>st'^) - 2r(tV + 18str - Qar'^)), 

= — ^-^(5t^ + 3(5^-^2) + 2tV2-3r^ + 6t(s2-a2 + 2(7r) 

-2<jr{e - 3r2) - {s^ - a^){be - 3r2)), 

= ^-^(-6tCT + tV + 6sr+ 18str + 3(s2 + a2)r-9CTr2 
2Ast^ 

+2r(-t2 + 3r2)), 
= — ^(3a-5tV + tV + 3ar2-3r3-6t(-(T + r)) . (3.11) 

With the denominator equation 1 — zst = 0, there are 5 equations in 5 variables. 
The system possesses only isolated solutions and both the Stringvacua package 
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and the NPHC method can solve this system. The CBB is 432, the 2HomBB, 
using {t, T, s, a} x {z} is 144, and the BKK root count is 100. There are 6 real 
solutions. 

Again, for a generic choice of the values for a, b and c, it becomes extremely dif- 
ficult for Stringvacua to solve the corresponding system unlike the NPHC method. 

Sys3-1: Two Moduli Example : We now take a four dimensional J\f = I 
supergravity theory with Kahler potential and superpotential 

K = -31og(Ti+fi)-31og(T2 + f2), 

W = -T^ - T1T2 - r| + 10(Ti + Ts) - 100. (3.12) 

Taking Tj = tj + iTj, for j = 1,2, the potential is 

^ = -^7^(ti+t?(20-14t2) + 20t^ + t^-60t2(100 + 2rir2 + r2') 

+t^(-500 + 280^2 - 37t^ + lOr^^ + 10rir2 + 7r|) 

+tl{7T^ + IOT1T2 + 10(-50 + r^)) - 2ti(-140t^ + 7tl 

+30(100 + T^ + 2T1T2) + 3t2(200 + r^ + 4rir2 + r^)) (3.13) 

+9(10000 + T^ + 2t^T2 - IOOt^ + 2rir| + r^ + t^{~100 + 3r|))) 

with the stationary points dictated by 

1^ = ^(-900ri-60tiri + 10t?ri-6tit2ri + 7t2ri + 18rf-60tir2 



+5tlT2 - Q0t2T2 - l2tit2T2 + 5t^r2 + 21t(t2 + 27rir^ + 9r| 



i^ = gMftl^"^^^!'"! + ^^i""! ~ ^°*2n - 12tit2ri + St^n + Qrt 
-900x2 + 7t?r2 - 60t2r2 - Qht2T2 + 10^^X2 + 27rf r2 + '. 
dV 1 / n^r\r\r\r\ 1 1 onnn-A 1 c;nn-A2 1 j.4 , 1 onnn+ 1 o/inn+ + ocn-A2 



ai 



dV 
dt2 



_j^ ^-270000 + 12000ti + SOOtf + t\ + 18000^2 + 2400tit2 - 280ttt2 
+1500^2 - hmtitl + 37^2^2 _ QQtl + 2Stitl - 2>tl 
+2700r2 + 120tirf - lOt^rf + 12tit2rf - 2ltlTl - 27Tf 
+2A0tiTiT2 - 10tlTiT2 + 360t2rir2 + 48tit2rir2 - 30tlTiT2 
-54TfT2 + 2700r| - 7t?r| + 180t2r| + 12tit2r| - 30tir| 
-81rfr| - 54rir| - 27r|), 
Y^^(-270000 + 18000ti + 1500t? - 60t? - 3tf 

+12000^2 + 2400tit2 - 560^2^2 + 28t?t2 + 500t| - 280tit^ + 37t?t^ + 4 
+2700r2 + 180tirf - 30t?r2 + 12tit2Tf - rtjr^ - 27Tf 
+360tirir2 - 30t?rir2 + 240^2^X2 + 48tit2rir2 
-10tlTiT2 - 5ATfT2 + 2700r| - 21tlT^ + 120t2r| 
+12tit2r| - lOt^rl - 81r2r| - 54rir| - 27r|). 

(3.14) 

We can easily clear out the denominators from each of the equations and 

add a constraint equation 1 — ztit2 = representing the fact that none of the 
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denominators is zero. This combined system contains only isolated solutions, but 
solving this system using the symbolic Grobner basis methods is a prohibitively 
difficult task |15| . 

However, using the NPHC method, we can solve this system in less than a 
minute on a standard desktop machine. The CBB of this system is 432, the 
2HomBB, using {ti, ri, ^2, ^2} x {z}, and the BKK root count are 144. In the 
end, there are 70 solutions, of which 6 are real. 

Sys3-2: Furthermore, the equations for F-terms are 



= 


= 300 - lOti -tl- 30t2 + tit2 + 3tl - 3rl - 3nT2 - 


- 3r|, 


= 


= 300 - 30ti + 3tl - iota + ^1^2 -tj- 3t^ - 3tiT2 - 


- 3tI 


= 


= -30ri + 2tiri + 3t2ri - 30r2 + tiTa + 6t2r2, 




= 


= -30ri + etiTi + t2Ti - 30t2 + 3tiT2 + 2t2r2 





The CBB and the BKK root count both are the same in this case, namely 16. Due 
to the lack of a good partition of the variables, the 2HomBB is actually larger, 
namely 24 for {^1,^2} x {ti,T2}, than the CBB for example. These equations 
have only 12 solutions, of which 4 are real. 

3.2 Examples from Heterotic and M-Theory Effective Po- 
tentials 

Having warmed up with some toy examples and seen that the NPHC method 
can extract information which is too difficult for Grobner basis methods or for 
ordinary numerical methods untuned for polynomial systems, let us see some 
"real" scenarios from phenomenology. 

Sys4: A Heterotic Example: Let us begin with a Kahler potential and 
superpotential taken from a heterotic compactification |1T| 



K = -31n(T + T)-31n(Z + Z), 

W = i{^ + ieT) + {e + ipT)Z + ^{fi + iqT)Z^ + l{p + irT)Z'^ , (3.15) 

2 o 

where ^, r, e, q, fi, p, p, e are parameters which satisfy the constraint ^r — eq + fip — 
pe = 0. We write T = t + ir and Z = z + i(, and take ^ = — 13,r = 0,e = 
— 4, g = 2,yU = 2,p = l,p = 5,e = —7 to ensure that the solution space is 
zero-dimensional. 
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The critical points satisfy 
= -ggj^ (164268 + 10152^2 _ 4932^^3 _^ 2268^^ 

-252*4^ + 360tz^ + Qlfyz^ + 101088^ - 11232^^ + 25380^2^ 
+1620z^C + 40824^2 _ i728tz^2 _^ 9720^2^2 _ 215*222^2 
+720tz3^2 _^ 2025z''^2 _^ 28836^3 + 3240^2^3 + 7452^^ 
+36*2^4 _^ 36otz^4 _^ 2025^2^4 _^ ig20^5 ^ 575^6 _ 8424z2t 
-216z^r - 4320z2^T - 1620z''C^ - 25272^2^^ _ 2160z'^^'^T 
-7776^3r - 3240z2^3j. _ 1944^4^. _ iq20^^T + lOSz^T^ 
+108022^2^.2 _^ 972^4j.2 _^ 1764(^2 _^ 27r2) 
-36(t2(7z2 - ^2) _ 3(^2 ^ 9^2)^2) ^ 84(-6t2e 
-6te(-24 + 25^2 - 12^ - 15^2) _ 27r(78 + 24^ + 6^2 + 5^3 

+^2(2 + 50 + 6CT) + 6(t2(llz2 + ^2) ^ 9(^2 ^ 3^2)^2)) 

-36(2tz(-78 + 6^2 + 10^3 _^ ^2(^4 _^ ^q^^J _ 3(5^4 

+^2(8 + 10^ + 20^2) _^ 3^(7g _^ 24^ + 6^2 + 5^^))T 

+2(t2(-3z2^ + ^3) + 3^(5^2 + 9e2)r2))), 

= -2304F?: (54756 + 1128^2 - 252^^ + 84*2^4 

-120tz^ _ 225^6 + 33696^ - 3744^^ + 2820^2^ - 180^^^ 
+13608^2 _ 576t2^2 _^ 1080^2^2 _ 72t2^2^2 

-225z^C^ + 9612^3 ^ 360^2^3 _^ 2484^4 _^ 35^2^4 
+120tz^'* + 225^2^4 _^ 54Q^5 _^ 225^*^ - 936z2r + 24z4r 
-480z2^r + 180z^Cr - 8424^2^^ _ 240^2^2 j. _ 2592^3^^ 
-360^2^32^ _ 648^42^ _ 540^5^. _ 12^4 j.2 

+ 12022^2^.2 _^ 324^4^.2 _^ 1764(^2 _^ 92^2^) 
+ 12(t2(-7z2 + 3^2) + (z2 + 27^2)2^2) 

+12(-2te(-78 + 6C^ + lOC^) - {5z'^ 

+z2(-8 - IOC - 20^^) - 9^(78 + 24^ + 6^2 + 5^^))T 

+2e(3t2(z2 _ ^2) _ (5^2 ^ 27^2)2-2)) 

+84(-6t2c + 6tz(8 + C(4 + 50) + 2(^2(11^2 + 3^2) 

+3(z2 + 9^2)2-2) _ 3T(72e + ^2(2 + 5^ 

+3(78 + 6^2 + 5^3 + 6er)))), 
= 3g4ij^(-3276 -240^2- 4^4 -1008^-290^2^-30^4^-720^^- 40^2^2 _ 354^; 

-60^2^3 _ 3g^4 _ 3Q^5 _^ ;^7g42. _^ I68z2r + 4z^T 

+504e2r + 40^2^22- + 3g^42- + 4(^2 ^ 9^2)^- 

+2(5^4 + ^2(8 + 010 + 20c - 20r)) + 3C(24C + 5^^ 
+3^2(2 - AT) - 6(-13 + UT)))), 
= 1152^3 3 (5616 - 936tz + 1410^2 

+90z^ + 4536C - 288tzS. + 1080^2^ - 72^2^2^ 
+120tz3^ + 2252^^ + 4806^^ + 540^2^2 
+1656^3 + 24*2^3 _^ i20tz^^ + 450^2^3 

+450^^ + 225^^ - 240z2r - 90z^T - 2808Cr - 240z'^^T 

-1296^'^T - 540^2^22^ _ 432^32^ _ 450^42^ 

+12022cr2 + 216^32^2 j_ i2^(t2 + 92^2) 

-42(-12tz - 4*2^ - 30tz^ + 72r + 15z2r + 36Cr + 45C^r 

-36CT2 + 2(t2 + 9T2)) + 2(-3(2tz(5z2 + 6^ + 15^2) 

+ (-144^ - ^2(10 + 400 - 3(78 + 18^2 j_ 20^3^)t] 

+6(3t2(z2 _ ^2) _ (5^2 ^ 27^2)2^2)))^ 
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(3.16) 



We can clear out the denominators from all the equations and add the additional 
constraint equation 1 — y zt = with additional variable y, to get a combined 
system which only possess of isolated solutions. The system is prohibitively 
difficult to solve using the Grobner basis techniques. However, the NPHC method 
solves it in less than a minute. The CBB is 2700, the 2HomBB using {t, r, z, (} x 
{y} is 900, and the BKK root count is 340. There are a total of 62 solutions, of 
which only 6 are real. 

Sys5: A Model from M-Theory: Here, we take an example of M theory 
compactified on the coset tt/s tt/-. |43| . The coset has SU(3) structure. The 
corresponding Kahler and superpotential are 

K = -41og(-z([/-f/))-log(-z(Ti-fi)(T2-f2)(T3-f3)), 
W = ^(4f/(Ti + r2 + T3) + 2r2T3-Tir3-TiT2 + 200). 



Using Tj = Tj — itj, for j = 1, 2, 3, and U = y — ix, the potential is 
^ = 7;777T^^(40000 + ^3^i -400rir2-4t^rir2 + 4t2r| + r2r2-400rir3 

+800X2X3 + 2t^T2T3 - 4TiT^T3 + T^T^ - ^T^T^tI + ^T^tI - 24^2^3^^ 

+4^3^^ - 24ti(t2 + h)x^ + AtIx^ + StiT2X^ + At^x"^ + 8rir3X^ + 8x2X30;^ 
+4:tIx'^ + leOOriy - Ulriy + 1600r2i/ + 16t^r2y - 8x^x22/ - ^titIv 
+1600x31/ - Srlny + IQrjny - Snr^y + 16r2r|y + 16^2^2 + 15^2^2 
+32TiT2y^ + 16r|i/2 + 32rir3y' + 32T2T3y^ + IGr^y^ 
+tl{tl + tl + tI + 2r2r3 + tI + 4^^ - 8x2?/ - 8x31/ + \^y'^) 
+^2(4^3 + ^1 - 4ri(r3 + 2y) + 4(r| + x^ + ^r^y + 4y2)). 

There are 9 equations to be solved: 8 equations corresponding to the deriva- 
tives with respect to 8 variables plus one constraint equation 1 — 2(^1^2^33^) = 
with z being an additional variable. The CBB is 103,680, the 2HomBB, using 
{tj,rj,x,y} X {z\ is 20,736, and the BKK root count is 18,624. The NPHC 
method yields 516 solutions, of which only 12 are real. The computation based 
on the BKK root count takes around 10 minutes using HOM4PS2. 

Sys6: A SU{2)xSU{2)/Z2xZ2 Model: Compactifying type HA supergravity 
on SU{2) X SU (2) /Z2 x Z2 = 5^ x S^/Z2 x Z2 and restricting to the modes that 
are left-invariant under the action of SU{2) x SU{2) one obtains an A/" = 1 
supergravity in 4D, which was studied in [HHIB] . There are no D-terms so 
that the scalar potential V is determined by the Kahler potential K and the 
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superpotential W which are 

K = - log((ti + ti)(t2 + i2){t3 + k)) 

- log(16(Ari + N,){N2 + N2)iNs + N3)iN^ + N^)) 
and 

W = -ifQhhh + tlt2/f ^ + tlts/f + hhf?^ + {-lh + h+t2+ h)Ni 

+ {-ih + ti-t2- h)N2 + {ih + ti-t2 + h)N^ + {-ih-h-t2 + ts)N4. 

Here tj are the Kahler moduh [J — iB = tjF^^"^*), A'^^ are the complex structure 
moduh [e~'^ ImiVt) + iC^ = NkY^^~^^^), /q is the mass parameter (of massive 
type IIA i.e. the flux Fq), fP is the Fa-flux (F2 = /f V^^^)^) and h is the in- 
flux (the F'4-flux is not turned on and also included 06-planes). The Y are the 
elements of the cohomology of M = SU(2) x SU{2)/'L2 x ^2- The superscript 
(i— ) means that they are elements of H^{M, R) and the superscript +/— means 
that they are even/odd under the orientifold projection, i runs over the number 
of odd 2-forms in cohomology. The fluxes and fields have fixed transformations 
under the orientifold projection and have to be expanded either even or odd 
cohomology in order to survive the orientifold projection. By rescaling the fields 
and the overall scale of V we can set /o, /} 5/2 , /a to unity. Furthermore, the 
scalar potential simplifies substantially if we analytically solve for the four C3 
axions (i.e. Im{NK))- This leaves us with a model having 10 real fields and one 
parameter h. 

As was shown in P5|H6]. the scalar potential V for this model circumvents 
all known no-go theorems against dS solutions and actually allows for dS criti- 
cal points that are however unstable. The results in PS | I ^ where obtained by 

numerically minimizing e = y2 using Mathematica. This method how- 

ever does not give all possible solutions. Using HOM4PS2, we have now solved 
d^V = and found all critical points for any given value of h. There are 10 vari- 
ables and 10 equations, plus an additional constraint equation with additional 
variable to make sure that the denominator is non-zero, yielding a total of 11 
equations. The CBB for this system is 279,936,000, the 2HomBB, using the orig- 
inal variables in one group and the additional variable in the other, is 34,992,000, 
but the BKK bound is only 574,080. For a fixed value of h{= —5) we have found 
that there are 60 real solutions to these equations. By stability analysis of the 
solutions, we conclude that only one (the previously known) unstable dS critical 
point exists in this case. Furthermore, we also find (previously unknown) that up 
to permutation of the fields there are 4 additional AdS solutions. Thus, we have 
now completely classified the stationary points for this model using the numerical 
homotopy continuation method. 
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Model 


# Eqns 


CBB 


2HomBB 


BKK 


Real Sols 


Timing 


GB Technique 


Sysl 


3 


364 


182 


182 


6 


< 20s 


Yes 


Sys2 


5 


432 


144 


100 


6 


< 10s 


Yes 


Sys3-1 


5 


432 


144 


144 


6 


< 30s 


No 


Sys3-2 


4 


16 


24 


16 


4 


< 10s 


Yes 


Sys4 


5 


2700 


900 


340 


6 


< 10s 


No 


Sys5 


9 


103680 


20736 


18624 


12 


~590s 


No 


Sys6 


11 


279936000 


34992000 


574080 


60 


11 hrs 


No 



Table 1: The short names of the various models used in the main text are written 
in the first column. The subsequent columns are the number of equations, CBB 
bound, 2HomBB bound, BKK bound, number of real solutions, and the time 
taken by the NPHC method (either Bertini or HOM4PS2). The final column lists 
whether or not the Grobner Basis technique (via Stringvacua) could solve the 
corresponding system. 



4 Numerical Algebraic Geometry 

We have presented, in the above, extensive examples wherein the power of numer- 
ical algebraic geometry may be harnessed in various phenomenological contexts, 
especially in the identification of isolated extrema of potentials coming from het- 
erotic and M-theory scenarios. One might naturally wonder as to situations where 
the solution space of algebraic systems has positive dimensionality. Such cases 
arise naturally as a central subject in the study of supersymmetric gauge theories. 

Indeed, the vacuum moduli space, parameterized by the scalar components of 
multiplets, in generic gauge theories with supersymmetry is a continuous man- 
ifold, or, more strictly speaking, an affine algebraic variety. This fact is par- 
ticularly pronounced in string theory when the gauge theory manifests as the 
world-volume theory on a brane and the vacuum moduli space is some space 
of special holonomy such as Calabi-Yau or G2 manifolds: this is the AdS/CFT 
correspondence. 

The geometrical engineering of world-volume gauge theories given an affine 
(Calabi-Yau) geometry has been a vast subject over the last decade (cf. e.g. |47] 
for an introduction) and conversely, the geometric analysis of the vacuum of a 
given supersymmetric gauge theory (whether it comes from string theory or not) 
could uncover hidden phenomenological symmetries [31I1S1I1S]- The algebraic ge- 
ometry of the vacuum could often be extremely complicated since the coordinates 
parameterizing it are the gauge invariant operators, subject to superpotential con- 
straints. Though in AdS/CFT, the geometries are affine Calabi-Yau threefolds 
by construction, reverse engineering [50] can lead to unlimited possibilities of 
Kahler manifolds. In super-symmetric QCD with gauge group SU{N^ with Nf 
fundamental flavours, for example, the vacuum is an affine variety of dimension 



as high as 2N^Nf - {N^ - 1) when Nj > N^. 

Furthermore, the vacuum moduh space could be composed of (possibly in- 
tersecting) unions of components - or branches - at various dimensions all the 
way from to that of the top (or coherent) component. This phenomenon is 
especially marked in the investigation of the so-called "master space", which is a 
join of mesonic and baryonic branches [5 H I52 ] . 

The importance of studying the algebraic geometric structure of the vacuum, 
as can be seen from the above-mentioned, should not be undermined by the tech- 
nical difficulties which we are inevitably lead to encounter. The usual method of 
attack is to recast the algebraic equations describing the vacuum as a polynomial 
ideal, transform into standard Grobner basis, and then perform primary decom- 
position to extract the irreducible branches. As is by now well recognized, this 
is a very expensive computation, given especially the double-exponential running 
time of computing Grobner bases (cf. [6]). In this section, we will study examples 
where if one wishes to know simply some crucial but preliminary information such 
as the number, dimension and degree of the components, one could bypass the 
prohibitive step of primary and irreducible decomposition and turn to the virtues 
of numerical algebraic geometry. 

For positive dimensional varieties, since there are infinitely many solutions, 
one needs a proper representation of the solutions. The important first question 
to be asked is how we need to represent the solutions. For a 0-dimensional variety, 
the solutions are just a finite set of points, so they can be represented by complex 
numbers. For positive dimensional varieties the situation is more involved. The 
solutions in this case form, for example, curves, surfaces, or hypersurfaces. A way 
to represent the solutions is to compute a parameterization of the varieties, which 
can be accomplished using a Grobner basis. However, this is computationally very 
expensive. The reader can find a nice discussion about this in [7]. 

Numerical Algebraic Geometry (NAG), on the other hand, cleverly uses an- 
other approach in which the solutions are represented as Witness Sets. We start 
with the fact that the number of points at which an irreducible component, say of 
dimension c? in C", of a system of polynomial equations intersects a random linear 
space of dimension c = n — d is equal to the degree of the component. For exam- 
ple, a cubic curve in C'^ intersects with a two-dimensional random hyperplane in 
three points. These intersection points are called Witness points. A random c- 
dimensional linear space means the solution set of d = n — c linear equations with 
random coefficients. Here, each linear equation is of the form Yl^=i ^j^j — C*o 
where Ck is a random complex number. After computing the Witness points, 
one can slide this random linear space around to obtain as many points on the 
component as needed. This yields a parameterization of the component, called 
the Witness sets representation. 

Now, it should be clear how the NPHC method plays an important role here. 
Algebraically, we add d linear equations to the system of equations and then 
compute the isolated solutions of the combined system using the NPHC method. 
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The number of isolated solutions gives the number of Witness points for dimension 
d, i.e., the degree of the union of components of dimension d. One can do this 
for all d between and n to compute the Witness points for all dimensions. We 
note that, in practice, a sequence of homotopies, called a cascade |53 | l5l ] . is used 
to compute the Witness points for all dimensions. 

Given the Witness points for the components of dimension d, one needs to 
partition these points into Witness points for the irreducible components. This is 
done by the so-called monodromy algorithm. The basic approach is to consider 
the Witness points as the linear space is moved around and then back to the 
original linear space. In this move, some of the Witness points return to the 
same point. However, if a Witness point returns to a different point, these two 
points must be on the same component and thus are grouped together. This 
process is repeated until each group is verified to be the set of Witness points for 
an irreducible component. Since the full technical description of this verification 
step of monodromy is beyond the scope of this article, we refer the reader to [29] 
for more details. 

A simple example to consider is that of the affine twisted cubic curve which 
is given by the equations: 

x^ -y = 0,x^ -z = 0. (4.17) 

Running Bertini, for example, on these equations yields that the irreducible de- 
composition of the variety contains one component of dimension 1 and degree 3, 
as expected. It also gives a random hyperplane that was used to slice the variety 
so that we can then get as many points as we like on this positive dimensional 
component. 

4.1 The Master Space 

First we apply NAG to a model where the results are at least partially known, so 
that we can have a bench-marking to our NAG set up. In [SUES], we introduced 
the concept of the "master space" of supersymmetric gauge theories, which con- 
trols the combined mesonic and baryonic branches of the vacuum moduli space. 
Computationally, this is the space of F-flatness, that is, the Jacobian ideal of 
the superpotential and can be construed as a baryonic fibration over the mesonic 
moduli space. In string theory, in the situation of a single D3-brane probing an 
affine toric Calabi-Yau threefold, the resulting world-volume physics is a U{1)^ 
quiver gauge theory; here, the master space is a toric variety of dimension g + 2, 
whose Kahler quotient by the g — 1 independent t/(l)-actions is the Calabi-Yau 
threefold. 

Though in the toric case, polytope and plethystic techniques can simplify 
the computation of the master space, direct attack thereon is still needed in 
general. Once again, for large number of fields, recasting the Jacobian ideal into 
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standard bases may become prohibitively difficult and numerical methods can 
become a more natural choice in quickly extracting the dimension and irreducible 
components of the master space. 

We will exemplify with A/" = 1 gauge theories for a D3-brane on the Abelian 
orbifold C^/Z^ x Z„. This is a toric case which had been approached by both 
direct and combinatorial methods for some low values of [m, n) in j51], and hence 
provide a good point of reference. Another advantage we have here is that the 
Jacobian ideal is always "square", in the sense that the number of variables is 
equation to the number of vanishing equations. To summarize, the gauge theory 
of interest is a quiver theory with mn nodes, a total of 3?7in bifundamental fields 
{Xjj, Yij, ^i.j} from node i to node j (with {i,j) defined modulo {k,m) and 
with superpotential 

k—l m—1 

W = 2_^ 2_^ ^*j^«+i,i-^i+ij+i ~ ^ij^ij+iZi+ij+i ■ (4-18) 

In Table |21 we provide the complete set of results for the systems up to 
{m,n) = (3,3). In jM], the top dimensional components of these systems were 
computed using toric variety methods. Our results agree with the results in this 
reference. In addition, we can now compute the complete irreducible decomposi- 
tion of all these systems hence already for this simple set of systems we have got 
new and complete results. 

For example, with {m,n) = (2,2), we arrive at a master space which is of 
dimension 16, in agreement with [51]. We further find that the variety reduces 
into 4 components: a dimension 6 piece of degree 14 and 3 linear pieces (i.e., 
degree 1) of dimension 4. 

One can also go further than these systems. For example, (4,1): 1 component 
of dimension 6 and degree 8; 33 components of dimension 4 and degree 1, (5,1): 
1 component of dimension 7 and degree 16; 131 components of dimension 5 and 
degree 1, and (6,1): 1 component of dimension 8 and degree 32; 473 components 
of dimension 6 and degree 1. For bigger systems such as (4,2), the decomposition 
space is more interesting, i.e., it has 1 component of dimension 10 and degree 
584, 9 components of degree 1 and dimension 8, 8 components of degree 3 and 
dimension 8, 8 components of degree 4 and dimension 8, 24 components of degree 
9 and dimension 8, and 8 components of degree 2 and dimension 6. 

For (5,2), there is 1 component of degree 3632 and dimension 12, 17 component 
of degree 1 and dimension 10, 10 components of degree 3 and dimension 10, 
40 components of degree 4 and dimension 10, 40 components of degree 9 and 
dimension 10, 40 components of degree 16 and dimension 10, 80 components of 
degree 27 and dimension 10, 130 components of degree 2 and dimension 8, 20 
components of degree 6 and dimension 8, and 20 components of degree 8 and 
dimension 8. We can go on for bigger and bigger systems, however, we would 
rather stop here and move towards a more difficult system next. 
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m\n 


1 


2 


3 


1 


NA 


(1;4;2 6), 
(1;2;1|6) 


(1;5;4|9), 
(7; 3; 19) 


2 


(1;4;2 6), 
(1;2;1|6) 


(1; 6; 14 12) 
(3;4;1|12) 


(1;8;92|18), 
(5; 6; 118) 
(6; 6; 3 18) 


3 


(1;5;49), 
(7;3;1|9) 


(1;8;92|18) 
(5; 6; 118), 
(6;6;3|18) 


(1; 11; 1620 27), 
(6;9;1|27), 
(27;9;2|27), 
(36;9;7|27), 
(27; 7; 1 27) 



Table 2: The master space F^ for C^/Z^ x Z„ as explicit varieties for several 
values of m and n. {p; q',r\d) means p components of dimension q and degree r 
contained in d-dimensional complex space. 

4.2 Supersymmetric Quantum Chromodynamics 

We now move to a more difficult problem of the familiar example of pure sQCD 
with gauge group SU{Nc), Nj number of fundamental flavours of quarks Q and 
anti-quarks Q, and no superpotential, so that the matter content is summarized 

as: 





GAUGE SYMMETRY 
SU(Nc) 


SU{Nf)L 


SU{Nf)n 


GLOBAL SYMMETRY 
Uil)B 


U{1)r 


U{1)q 


Uii)a 


Ql 


□ 


D 

1 


1 
□ 


1 
-1 


Nf 
Nf-N, 

■iv, 


1 






1 



(4.19) 

The generators of the gauge invariant operators consists of the mesons and 
baryons: 

M'j = QiQ'^ (mesons) ; 



ll...tNc ^ll • • • ^iNc "1 



ai...ajVc 



fi 



■CLNc 



(baryons) ; 
(antibaryons) . 



(4.20) 



It is a standard fact that the dimension of the (classical) vacuum moduli space 
Ai(Nf,Nc) of ths theory is 



dimA^ 



iNf,N,) 



^h 



Nf<N,- 
2N,Nf - {Nl - 1) , Nf>N,. 



(4.21) 



More speciflcally, one can explicitly obtain the algebraic variety [l9] and prove 
they are, in fact, all affine Calabi-Yau varieties (the reader is referred to Table 
2 in |in]); for example, Ai{Nf,Nc) is simply the affine space C f when Nf < Nc, 
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Nf\N, 


1 


2 


3 


4 


1 


(1;2;23) 


- 


- 


- 


2 


(1;4;6 8) 


(1;5;2 6) 


- 


- 


3 


(1;6;20 15) 


(1;9;14 15) 


(1;10;3 11) 


- 


4 


(1; 8; 70 24) 


(1;13;132|28) 


(1;16;115|24) 


(1; 17; 4 18) 



Table 3: The irreducible decomposition of the moduli space of sQCD the- 
ories with Nf fermions and Nc bosons (after elimination from Macaulay 2). 
Here, Nf > Nc, and the irreducible decomposition is presented in the form: 
(no. components; dimension; degree | total dim). 

and, for A^^- = N^ it is a complete intersection, with the Hilbert series given as 



gNf=N^ 



it) 



1-t 



2Nc 



{1 - t^)^^i - t^^y 



(4.22) 



In general, however, the vacuum is quite involved. Primary decomposition to 
find the irreducible branches of the moduli space, as was experimented in |49| . 
is prohibitively difficult, even for small values of Nf and N^. It was conjectured, 
by working over the coefficient fields of the rationals or number fields of finite 
characteristic, that Ai{Nf,Nc) is actually irreducible. Now, we can properly check 
the components by working over the complex numbers, as one should. For exam- 
ple, take {Nf,Nc) = (3,3), we readily find that there is only one component, of 
complex dimension 10 and degree 3. For [Nf, Nc) = (4, 3), there is again only one 
component of dimension 16 and degree 115. A full study of the vacuum structure 
requires algebraic elimination, a subject which we will address systematically in 
a forthcoming publication ^^ . 

For now, we can extract information about the full mesonic moduli space by 
allowing Grobner techniques to perform the elimination, such as using Macaulay2, 
and then use our numerical methods to perform the hard step of primary decom- 
position. We tabulate some results in Table [31 In ^\, only the top components 
of these systems using Macaulay2 were obtained and conjectured that the results 
were irreducible. In our case, we can now get the complete irreducible decom- 
position, not only the top dimension, of these systems. Hence our results yield 
that the conjecture was correct. We could have easily been able to go beyond the 
(4,4) systems, but the limiting factor was Macaulay2 itself which failed to obtain 
the eliminated ideals for the bigger systems for us. 

Indeed, one could combine our present course of study with that of the pre- 
vious subsection. Indeed, the standard Yukawa term to sQCD for Nf > Nc, in 
the above notation, is 



Nf Nf~Nc 
i,j=l a,b=l 



(4.23) 
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Nf\N, 


1 


2 


3 


4 


1 


(1;0;12) 


- 


- 


- 


2 


(1;0;14) 


(1;5;4|8) 


- 


- 


3 


(1;0;16) 


(1;7;6 12) 


(2;11;18|18) 


- 


4 


(1;0;18) 


(1;9;8 16) 


(2;14;32|24) 


(2;19;88|32), (1;20;320|32) 



Table 4: The master space of sQCD theories with Nj fermions and Nc bosons. 
Here, Nf > Nc, and the irreducible decomposition is presented in the form: 
(no. components; dimension; degree | total dim). 

Here, Mjj is the meson condensate of the quark-anti-quark. This superpotential 
can be considered as being generated by Seiberg duality from one with pure 
matter content. For example, at {Nf,Nc) = (3,2), we find the master space to 
be of dimension 9, with 3 components, one at degree 1 and two at degree 7. The 
complete results are shown in Table [Hand we see non-trivial primary components 
to the master space. 



4.3 Instanton Moduli Spaces 

The study of Yang-Mill instantons is, undoubtedly, another important subject 
in the investigation of gauge theories. Ever since their discovery by Belavin, 
Polyakov, Schwartz and Tyuplin [HS] as well as the construction by Atiyah, Drin- 
feld, Hitchin and Manin (ADHM) on self-dual solutions j^^ in the 1970's, the 
parameter, or moduli, space of these instanton solution has been of great atten- 
tion to physicists and mathematicians alike. The geometry of the moduli space 
can be quite involved. Indeed, whereas the ADHM construction gives the moduli 
space for the classical Lie groups, that for the exceptional ones still remain a 
mystery. 

With the embedding of the construction into string theory [58] by Witten, 
Douglas and Moore, one could realize many cases as the vacuum moduli space of 
the supersymmetric (quiver) theories which we have introduced above. Recently, 
the algebraic geometry of the one-instanton moduli space was analyzed along the 
line of our present thought in [52] , calculating, specifically, the Hilbert series of 
the space. To clarify notation, we shall let k G-instantons signify instantons of 
(classical) gauge group G and with winding number k. Let us take Figure 7 of the 
said paper as the illustrative starting point. Here, the moduli space of k SU{N) 
instantons is given as the A/" = 1 vacuum moduli space of the "flower quiver", 
with the field content: 0(*=^'^) and $ charged as k x k matrix fields under U{k), 
X21 charged as N x k and X12 as k x N matrix fields respectively, all obeying 
the superpotential: 



ly = X21 • $ • X12 + e„fl0^"^ ■ $ ■ <^^"^ . 



(4.24) 
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Since the case oi k = 1 was considered in great detail in ^B\, let us move on 
the much more complicated and unsolved problem of, say, k = 2. We will focus 
on the Higgs branch by setting $ to have zero vacuum expectation value. The 
input data to the moduli space is as follows. 

Gauge invariants: 

^*=l / ji,J2=l,-,N 

2 



1^ ^12 (PhA2^2 



F- Terms: 



■12 V^Ji,J2^'-21 



= %elds k=o [Tr X21 ■ $ ■ X12 - Tr $[0(i), 0(2)]] . 

(4.25) 
Let us be more specific and take, for example, at A^ = 2. First, the master space 
is easy to determine. It is a dimension 12 variety of degree 16, defined by 4 
quadrics in C^^. The Hilbert series is 



9{t;^\k,N)=i2,2)) = TTZlk PE~'[9m = m-Ae . (4.26) 



{1-ty 

From the termination of the plethystic logarithm [601, we see that it is, as also 
suggested by the dimension, a complete intersection. 

This is an ideal with 21 generators in 33 variables: 16 variables in the X 
and fields together with 17 auxiliary variables. It can be easily shown that 
this ideal has a unique irreducible component of dimension 12 in C^^, but com- 
puting its degree and the actual moduli space via elimination are already quite 
overwhelming for standard computer algebra packages. Bertini, in just under 24 
hours using 200 processors in parallel, found that the degree of this irreducible 
component of dimension 12 is 20364. We will report on computing the moduli 
space via numerical elimination in a forthcoming paper. 

5 A Few Words About NAG Packages 

Although it is not at all our intention to compare the performances of different 
packages in this paper, we convey what the important computations are that 
different packages can do more efficiently than others. There are so far three 
independently packages for NAG: Bertini, HOM4PS2 and PHCPack. 

Bertini is a general purpose package for NAG, i.e., it can find all the stationary 
points of a 0-dimensional ideal and also find irreducible decomposition of positive 
dimensional ideals efficiently. It also has many other Numerical Algebraic Geom- 
etry implementations such as a membership test, finding multiplicities of isolated 
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singular solutions, facility to specify user-defined honiotopy, etc. In particular, 
Bertini automatically constructs Total Degree and 2- Homogeneous Homotopies. 
It utilizes adaptive multiprecision, and a parallel version of Bertini is also publicly 
available. 

HOM4PS2 is regarded as having the fastest path-tracker |36]. It also has both 
Total Degree and Polyhedral Homotopies. HOM4PS2 uses MixedVol-2.0 which 
is a very efficient package to compute the BKK bound. On the other hand, the 
applications of the current version of HOM4PS2 are limited: for the 0-dimensional 
case, it only takes one of the two homotopies as its only input parameter along 
with the equations themselves and gives out all the isolated real and complex 
solutions. For the positive dimensional systems, it also takes in only two of the 
homotopies as its input parameter and the dimension of which one wants to 
get the irreducible decomposition of along with the equations, and returns the 
irreducible decomposition of that particular dimension. The parallel version of 
HOM4PS2 is under development and will be soon publicly available |61] . 

PHCPack is a general purpose solver like Bertini. It does quite a few things, 
such as computing many different upper bounds on the number of solutions in 
addition to the CBB, 2HomBB, and BKK root counts, it also can preprocess 
the system before running the path-tracker on it, one can supply a user-defined 
homotopy as well, one can choose different numerical precision and parameters 
manually if one wants to, etc. It runs interactively and one can chose different 
options in its menu-based interface. The parallel version of PHCpack is also 
publicly available. 

All these packages have their own advantages. It is a goal of the authors of 
the current paper to write a separate article on demonstrating various features of 
these packages specifically to the string phenomenologists and particle physicists, 
and write a Mathematica interface for these packages in the same fashion as 
Stringvacua. 

6 Conclusion and Prospects 

With the advances in computer algebra and algorithmic geometry, we have wit- 
nessed an increasing trend over the last decade wherein such a paradigm has 
been gaining importance in theoretical physics, especially in gauge and string 
theories. Not only can hitherto unthinkable problems such as stabilizing highly 
complicated potentials or scanning through huge classes of string vacua become 
feasible, but also new geometrical quantities which characterize gauge theories 
such as Hilbert series or topological invariants of vacuum moduli spaces can be 
calculated to elucidate the physics. 

One great hurdle to this approach of computational algebraic geometry is that 
the central method involved is finding the Grobner basis of polynomial ideals, a 
highly non-parallelizable and exponential growth algorithm. In this paper, we 
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have demonstrated many instances where numerical algebraic geometry, espe- 
cially using homotopy continuation methods, can supplement this short-coming 
when addressing certain quantities. Though by no means a replacement for 
the powers of the Grobner techniques due to its symbolic nature, the numeri- 
cal methods presented in this paper are very useful and often computationally 
less expensive. In addition, the numerical methods are parallelizable making it a 
spectacularly robust tool to address the string phenomenology related problems. 

Two classes of problems, arising from a myriad of physical situations such as 
effective potentials in string and M-theory compactifications and supersymmet- 
ric vacua of gauge theories, are particularly amenable to our method: finding 
critical points of zero-dimensional ideals and irreducible decomposition of higher- 
dimensional polynomial ideals. The former is a typical problem in vacuum sta- 
bilization problems and the latter, in extracting branches of the moduli space of 
vacua. We have shown many concrete examples where these are beyond current 
computational powers using Grobner bases but are fairly quickly done, in a highly 
parallelizable fashion, using numerical homotopy. 

We are clearly only touching the surface of a reservoir of great utility. The 
two classes of problems we have focused on in this paper, already of diverse appli- 
cations, are only two of the many quantities and techniques which have recently 
emerged to be crucial in the study of gauge and string theories. For example, the 
computation of the Hilbert series of a polynomial ideal, a problem seemingly of 
interest only to pure algebraic geometry, has turned out to enumerate supersym- 
nietric BPS spectra of operators |6U| . Schubert calculus is another fascinating 
area coming from enumerative geometry where the NAG method can be directly 
and efficiently applied. A new homotopy meant to deal with the Schubert calcu- 
lus problems has already been worked out. Currently, methods for numerically 
obtaining the Hilbert series are being developed. As another example, systematic 
elimination of variables in the so-called syzygy problem has turned to be impor- 
tant in the study of moduli space of vacua in gauge theories and give explicit 
defining equations of the vacuum manifold j3]. This, too, is being developed. 
Indeed, there are several international collaborations which are engaged in the 
combing over the plethora of string vacua using computer search. Most of these 
ultimately reduce to manipulation of polynomial ideals and are thus limited in 
parallelizability and running time by the Grobner paradigm. The possibilities of 
numerically computing such sophisticated quantities as Hilbert series and ranks 
of (co-)kernels of polynomial maps should offer a new outlook. 
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A 2-Hoinogeneous Homotopy 

As seen above, the variables in a polynomial system can often be partitioned 
into two groups. The key is a refined notion of the degree of a polynomial. To 
illustrate, consider the polynomial f{y, z) = 1 + 2y — 5z + 3y^z. Clearly, the 
degree of / is three. However, y has a maximum degree of 2 and z only appear 
linearly in /. That is, if a and /3 are random constants and 

9{y) = f{y, /3) and h{z) = f{a, z), (A.27) 

the degree of (7 is 2 while the degree of /i is 1. We say the bidegree of / is (2, 1). 

The notion of the bidegree of f{y, z) naturally extends to the case when y and 
z are sets of variables. That is, the bidegree of f{y, z) is (a, b) if the degg = a 
and deg h = b where g and h are defined by (1A.27P with a and /3 random vectors 
of the appropriate size. 

For a system P{x) = P{y, z) of m polynomials in m complex variables, where 
y is k complex variables and z is m — k complex variables, let (oj, bi) be the 
bidegree of Pi{y,z). When each Pi is nonconstant, the polynomial 

m 

R{s,t) = Y[{aiS + kt) (A.28) 

is homogeneous of degree m. The 2-Homogeneous Bezout Theorem (see (25] for 
more details) states that the number of isolated solutions of P{x) = in C™ is 
bounded above by the coefficient of the s^f^'^ term of R. This coefficient is the 
2-Homogeneous Bezout bound (2HomBB). 
For example, consider the system 

h{y,z) = 3y-2z + 4z^ = 0, 

h{y,z) = l-yz = Q, (A.29) 
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where y and z are complex variables. The CBB yields that this system has at most 
3-2 = 6 isolated solutions in C^. The bidegree of /i(?/, z) and /2(y, z) is (1, 3) and 
(1,1), respectively, with the polynomial i?(s, t) = (s + 3t)(s + t) = s^+4st + 3t^. 
Since the coefficient of st in R is 4, the 2HomBB yields that (1A.29|) has at most 
4 isolated solutions in C^. 

To construct a 2-Homogeneous Homotopy, we need to describe how to con- 
struct the start system Q{x) = Q{y,z). The polynomials Pi{y,z) and Qi{y,z) 
must have the same bidegree, namely (cj, 6j). One way to construct such a poly- 
nomial is to take 



Q.{y,z) = UUM ■ \\M,,{z) (A.30) 




where Lij{y) and Mij{z) are random linear polynomials. The number of solutions 
of Q{x) = Q{y,z) = is the 2HomBB of P{x) = P{y,z), and the solutions 
themselves can be computed using linear algebra. 

B Polyhedral Homotopy 

As mentioned above, polynomial equations arising in real-life problems are sparse 
in terms of the number of monomials which appear. In general, CBB becomes 
an upper bound on the number of solutions for such cases. We are looking for 
a tighter bound that takes the sparsity of the system into account. There has 
been a huge amount of work done on related issues using resultants and algebro- 
geometric methods, but the most important result for us is Bernstein's theorem. 
In order to state it clearly, we remind the reader of some standard notions of 
Laurent Polynomials, Newton Polytopes, and Mixed Volume. 

A Laurent polynomial allows negative exponents for the monomials, so no 
variable is allowed to be zero. Hence, multiplication of any Laurent polynomial 
by a monomial does not change the root count in (C*)"* = (C/jO})™". Formally, 
let Si C Z"* be a set of vectors whose elements are the exponents of the monomials 
of the ith polynomial. Si is called the support of the ith polynomial. Then a 
polynomial (say, ith polynomial) of the form fi{x) = Xlaes ^j^qX" is called a 
Laurent polynomial. Here, q^q, G C are the coefficients of the monomial x" with 

X e (C*)"*. 

Next, a set of points is called a convex set if for every pair of points within 
the set (or more formally, the mathematical object made by the set), every point 
on the straight line segment that joins them is also within the set. The convex 
hull of a set X is the minimal convex set containing X. We note that the convex 
hull of support Si of a polynomial, say Qi = conv(S'j), is called the Newton 
polytope oi fi{x). 
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For example, consider a two-variable system 

fi{x,y) = 1 + ax + bx'^y'^ = 0, 

f2{x,y) = 1 + ex + dy + exy'^ = 0, (B.31) 

where x and y are complex variables and a, b, c, d, e are complex coefficients. 
The CBB of this system is 4 • 3 = 12, i.e., there can be a maximum of 12 
isolated solutions for this system in C^. The 2HomBB, using {x} x {y} is 6. 
Now, the supports of these equations are ^i = {(0, 0), (1, 0), (2, 2)} and 5*2 = 
{(0, 0), (1, 0), (0, 1), (1, 2)} respectively. The Newton polytope for fi{x, y) is Qi = 
conv(5i) = {(0,0), (1,0), (1,1), (2, 2)} and for h{x,y) it is Q^ = conv(52) = 
{(0,0), (1,0), (0,1), (1,1), (1,2)}. 

A Minkowski sum of any two convex sets is defined as 

Qi + Q2 = {gi + g2 : gi G Qu q2 e Qi}- (B.32) 

The Minkowski sum of two Newton polytopes corresponds to multiplying the 
corresponding polynomials algebraically. The m-dimensional volume of a simplex 
having vertices vq, Vi, . . . , Vm, is 

Yolmiconv{vo, . . . ,Vm)) = — |det[t;i -vo,...,Vm -Vo]\. (B.33) 

From there on, one can show that the ?7i-dimensional volume Vo\m{)^iQi + ■ ■ ■ + 
XmQm), where > Aj G M, is a homogeneous polynomial of degree m in variables 
Aj. The mixed volume of convex polytopes Qi, . . . , Qm, denoted M{Qi, . . . , Qm), 
is defined as the coefficient of Ai ■ ■ ■ A^ in Vol„^(AiQi + ■ ■ ■ + XmQm)- It can be 
shown that 

m 

M{Q,, ...,Qm) = Y.^-ir-' VoU( Y, Qi)^ (B.34) 

where the inner sum is a Minkowski sum of polytopes and fi™ are the combina- 
tions of m-objects (i.e., m-dimensional geometrical objects made of m-simplices) 
taken i at a time. Moreover, the mixed volume is always an integer for a system 
of Laurent polynomials. For the case of two polynomials in two variables, 

M(Qi,Q2)= Vol2(Qi + Q2)- Vol2(Qi)- Vol2(Q2). (B.35) 

For the system in flB.3ip . 



voi2(gi) = 1, 

Vol2((52) = area of parallelogram made by {(0, 0), (1, 0), (0, 1), (1, 1)} 
-|- area of triangle made by {(1, 1), (1, 0), (1, 2)} 

= ^+2 = 2' 
13 

voi2(gi + g2) = —. (B.36) 
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Thus, the mixed volume for this system is 4. This is important because the 
Bernstein theorem and Bernstein-Khovanskii-Kushnirenko (BKK) theorem fS7\ ~ 
[39] says that for generic coefficients, the number of isolated solutions in (C*)™' of 
a Laurent system is exactly equal to the mixed volume of this system counting 
with multiplicity. For any particular set of coefficients, this is an upper bound. 
This result is very interesting since one can get a generically sharp bound on the 
number of solutions in (C*)™ of a polynomial system by knowing the exponent 
vectors of monomials. For the system in flB.3ip . there can be a maximum of 4 
isolated solutions in (C*)^, which is in contrast to the CBB and 2HomBB bound 
of 12 and 6, respectively, in C^. 

To yield the BKK count for C™, one considers an extension of the mixed 
volume, called the stable mixed volume\j, which ensures that we have all necessary 
solutions in C"^ ISSHHl]- A discussion on the stable mixed volume is beyond the 
scope of this article, although its calculation is similar to that of the mixed volume 
with some formal complications. However, it should be noted that a highly 
sophisticated implementation of an algorithm to calculate the mixed volume of a 
given system is MixedVol [33] which is transplanted in PHCpack and MixedVol-2.0 
which is transplanted in HOM4PS2. 

Now, after calculating the stable mixed volume of the original system P{x) = 
0, we require a start system, Q{x) = 0, such that Q{x) has the same stable 
mixed volume, where of course solutions of Q{x) = should be known or can be 
obtained easily. This homotopy method is called Polyhedral Homotopy. 



^Because of the highly technical nature of the stable mixed volume, the bound in C™ which 
has been commonly used and implemented in the community is the bound given in [S2]. This 
bound is quite easy to state: Add a constant term to polynomials in the system which do not 
have constant term, and the mixed volume of the resulting augmented system serves as a bound 
in C"' for the original system. The "stable mixed volume" is a little bit more general then this 
"augmented mixed volume". We would like to thank T.Y. Li for clarifying this point. 
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